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Abstract — Establishing bounds on the accuracy achievable by 
localization techniques represents a fundamental technical issue. 
Bounds on localization accuracy have been derived for cases in 
which the position of an agent is estimated on the basis of a set of 
observations and, possibly, of some a priori information related 
to them (e.g., information about anchor positions and properties 
of the communication channel). In this manuscript new bounds 
are derived under the assumption that the localization system is 
map-aware, i.e., it can benefit not only from the availability of 
observations, but also from the a priori knowledge provided by 
the map of the environment where it operates. Our results show 
that: a) map-aware estimation accuracy can be related to some 
features of the map (e.g., its shape and area) even though, in 
general, the relation is complicated; b) maps are really useful in 
the presence of some combination of low signal-to-noise ratios 
and specific geometrical features of the map (e.g., the size of 
obstructions); c) in most cases, there is no need of refined maps 
since additional details do not improve estimation accuracy. 

Index Terms — Localization, Cramer-Rao bound, Ziv-Zikai 
bound, Weiss-Weinstein Bound, A Priori Information, Map. 



I. Introduction 

Recently, localization systems have found widespread applic- 
ation, since they allow to develop a number of new services in 
both outdoor and indoor environments (TJ, ||2). Conventional 
localization systems acquire positional information from a set 
of observations; these are usually extracted from noisy wireless 
signals propagating in harsh environments. Unavoidably, this 
limits the accuracy that can be achieved by such systems. 
Localization accuracy can certainly benefit from the availa- 
bility of any form of a priori knowledge. In the technical 
literature, the sources of prior knowledge commonly exploited 
are represented by the positions of some specific nodes (called 
anchors) or by some characteristic of the communication 
channel; for instance in [3 j-[6| the impact of different factors 
(like non-line-of-sight, NLOS, propagation and network syn- 
chronization) on localization accuracy have been thoroughly 
analysed employing Cramer-Rao bounds (CRBs) and, when a 
priori knowledge is available, Bayesian Cramer-Rao bounds 
(BCRBs). In [7|-[10| multipath propagation in ultra wide 
band (UWB) systems and its effects on time-of-arrival (TOA) 
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estimation have been investigated by means of CRBs or 
BCRBs. In (TT)-(13) CRBs and BCRBs for the analysis of 
cooperative localization techniques have been derived. 
In recent times, however, some attention has been paid to the 
possibility of improving accuracy by endowing localization 
systems with map-awareness, i.e., with the knowledge of 
the map of the environment in which they operate. Specific 
examples of map-aware algorithms have been developed in 
I|l4|-|fl8|, which propose the adoption of non-linear filtering 
techniques (namely, particle filtering and extended Kalman fil- 
tering) to embed map-based a priori information in navigation 
systems. This work evidences that this kind of information 
plays an important role in improving localization and navi- 
gation accuracy; however, the impact of map-awareness on 
the performance limits of localization systems is still an open 
problem. In fact, as far as we know, the technical literature 
dealing with the fundamental limits of localization accuracy 
(9), pO) , p2) , fl3) has not considered this issue yet. 
In this paper, novel accuracy bounds for map-aware localiza- 
tion systems are developed and their applications to specific 
environments are analysed. Specifically, the BCRB, the exten- 
ded Zik-Zakai bound (EZZB) and the Weiss-Weinstein bound 
(WWB) for the above mentioned systems are derived. These 
bounds, which provide several novel insights into map-aware 
localization, have the following features: 

1) they are characterized by different tightness/analytical 
complexity; 

2) they can be evaluated for any map geometry and, in the 
cases of the BCRB and the EZZB, admit a closed form 
for rectangular maps; 

3) they allow to a) identify map features (e.g., shape, area, 
etc.) influencing localization accuracy and b) quantify 
such an influence. 

Then, such bounds are evaluated and compared for simple 
representative environments. Our results allow to assess the 
importance of map-awareness in localization systems in the 
presence of noisy observations. In particular, they evidence 
that: 

1) maps should be expected to play a significant role at 
low signal-to-noise ratios (SNRs) and in the presence of 
obstructions; 

2) in some cases simplified map models can be adopted for 
localization purposes, since map details have a negligible 
impact on estimation accuracy. 

The manuscript is organized as follows. In Section [IT] our 
reference scenario is illustrated and the map model is defined. 
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Performance bounds for map-aware localization are derived 
and evaluated in specific scenarios in Sections [III] and |IV] 
respectively. Finally, some conclusions are offered in Section 

m 

Notations: Matrices are denoted by upper-case bold letters, 
vectors by lower-case bold letters and scalar quantities by 
italic letters. The notation E r { } denotes expectation with 
respect to the random vector r; tr{-} denotes the trace of 
a square matrix; diag{-} denotes a square matrix with the 
arguments on its main diagonal and zeros elsewhere; is 
the element on the i-th row and Z-th column of its argument. 
Ijv denotes the N x N identity matrix; A >z B means that the 
matrix A — B is positive semi-definite. The probability density 
function (pdf) of the random vector R evaluated at the point 
r is denoted as /(r), whereas A/"(r;m, X) denotes the pdf 
of a Gaussian random vector having mean m and covariance 
matrix X, evaluated at the point r. I5 (x) denotes the so called 
indicator function for the set S (it is equal to 1 when x G S 
and zero otherwise). 

II. Reference Scenario and Map Modelling 

In the following we focus on the problem of localizing 
a single device, called agent, in a 2-D environment, i.e., of 
estimating its position p = [x,y] T , in the presence of the 
following information: a) a map representing the a priori 
knowledge about the environment in which the localization 
system operates; b) an observation vector z related to the true 
position p of the agent. 

A. Map Modelling 

In a Bayesian framework the map of a given environment 
can be modelled for localization purposes as a pdf f(p) of the 
random variable p, representing the agent position. Such a pdf 
is characterized by a 2-D support TZ having a finite size and 
consisting of the set of points in the environment not occupied 
by obstructions (e.g., walls, buildings, etc). In the absence of 
other prior knowledge, a natural choice for /(p) is a simple 
uniform model, i.e., 



/(p) = /(*,!/) 



l/A n pen 
elsewhere 



(1) 



where An denotes the area of TZ. This model is referred to as 
uniform map in the following and is preferable to other, more 
detailed, prior models since it can be adopted when only basic 
knowledge of the environment is availabl^] (e.g., the floor 
plan of a building in indoor environments or satellite photos 
in outdoor environments). An alternative statistical model for 
maps, better suited to mathematical analysis in the BCRB case, 
is a "smoothed" version of ([T}; this is obtained modifying the 
uniform pdf in a narrow area around the edges of TZ in order 
to introduce a smooth transition to zero along the boundaries 

'More detailed prior modelling requires additional information both for 
outdoor and indoor environments (e.g., the end use of rooms, the authoriz- 
ation levels to access them, furniture disposition, etc). Such information are 
often time-variant (e.g., the furniture may be moved), agent-specific (e.g., 
authorization levels and human habits may vary) and difficult to acquire. On 
the contrary, uniform maps only require the knowledge of building maps. 



of its support. Note that the adoption of a smoothed uniform 
model leads to analytical results similar to those found with 
the uniform model if the transition region is small with respect 
to the size of TZ. 

As it will become clearer in the following, for a generic 
map, a detailed description of the structure of its support TZ 
should be provided to ease the formulation of the accuracy 
bounds referring to the map itself. 

To begin, let us define, for each (x, y) <E R 2 , the sets 
7Z h (y) = {t\(t,y) € TZ} and TZ v {x) = {t\{x,t) € TZ} repres- 
enting the intersection between TZ and an horizontal (h) and 
vertical (v) line identified by the abscissa x and ordinate y, 
respectively. Then, let: 

• Nh(y) and N v (x) denote the number of connected com- 
ponents of TZh{y) and TZ v {x), respectively; 

• Af%(y) and M°(x) denote the set of the odd values 
{1, 3, 2N h {y) - 1} and {1, 3, 2N v (x) - 1}, re- 
spectively; 

• w n{y) an d h m (x), with n e A/^(y) and m e N°{x), 
denote the length of the connected components of TZh (y) 
and TZ v (x) respectively; 

• Af%(y) and Af^(x) denote the set of the even values 
{2A,- : ,2(N h (y)-l)} and {2,4,..., 2(N v (x)-l)}, 
respectively. 

• Aw n (y) and Ah m (x), with n <G Af%(y) and m <E N^{x), 
denote the length of the connected components of TZ\ (y) 
and TZ%(x) respectively (note that TZ^y) and TZ% (x) rep- 
resent the complementary sets of TZh{y) and TZ v (x) with 
respect to K); in other words, they denote the horizontal 
and vertical size of map obstructions, respectively; 

• X and y denote the projection of TZ on the x and y axis, 
respectively; 

The graphical meaning of these parameters for a specific 
uniform map is exemplified by Fig. [T] 

B. Observation Modelling 

In the following, we assume that the localization system 
is able to acquire a set of noisy observations, collected in 
a vector z and related to the agent position p according to 
a specific statistical model. Note that, on the one hand, the 
results illustrated in Sec. |III-A| hold for a general statistical 
model relating z to p; on the other hand, in Sec. |III-B| and 
|III-C| it is assumed that 



/(z|p)=AT(z;p,S) 



(2) 



for mathematical convenience, where z e R 2 and X = 
diag { fi 2 , <7 2 } (i.e., z is affected by uncorrected noise). It 
is important to point out that: 

1) The model |2|i is purposely abstract and does not ex- 
plicitly refer to a specific localization technique (e.g., 
TOA or received signal strength (RSS)) or to particular 
propagation conditions. Generally speaking, it is suitable 
to describe the position estimate z generated by a map- 
unaware and unbiased localization algorithm; in addi- 
tion, it leads to useful bounds which unveil the impact 
of map modelling (instead of that of observation model- 
ling) on estimation performance; readers interested in an 
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N v (x) = 1 
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Figure 1. Support region 1Z = IZi {JTZ2 of a non-convex uniform map (the notation introduced in Sec. |II-A] has been adopted). Note that 1Z (grey region) 
represents the area where the agent has non-zero probability to lie, so that everything else has to be considered as an "obstruction" (see fTJ). 



in-depth analysis on the impact of observation modelling 
can refer to [3|-[10|. 
2) In principle, the bounds derived in this manuscript can 
be extended to the case of correlated noise, accounting 
for the presence of an off-diagonal term <j xy ^ in S. 
However, the presence of noise correlation is neglected 
in the following derivations since: a) the magnitude 
of a xy in a real world system depends on several 
parameters (the type of measurements processed by 
the system, the multi-lateration or angulation technique 
employed, the properties of the propagation channel, etc) 
and cannot be easily assessed; b) our research work 
evidenced that the presence of noise correlation does 
not provide significant additional insights on the impact 
of map-awareness on performance bounds. 

III. Performance Bounds for Map-Aware Position 
Estimators 

In this Section various bounds about the accuracy of map- 
aware position estimation are derived and their implications 
are analysed. 

A. Bayesian Cramer-Rao Bounds 

Given an estimator p(z) of p based on the observation 
vector z, the BCRB establishes that the Bayesian mean square 
error (BMSE) matrix 

A(p(z))^E z , p {(p(z)-p) (p(z)-p) T } 

satisfies A(p(z)) t J 1 (H)}, Ipol, where J is the Bayesian 
Fisher information matrix (BFIM) associated to the statistical 
models employed by p(z). This result entails that, for 2-D 
localization: 



E 4 diag {E x ,E y } hB^ diag {B x , B y } 



where E x 4 [A (p(z))] M , E y 4 [A (p(z))] 2 2 , B x 4 
and By = [J _1 ] 2 2 - It is not difficult to show that 
the BFIM J can be put in the form (e.g., see fl9] p. 183, eq. 
(75)]) 



•z|p 



where 



and 



<z|p 



E* 



d_ 
dp 



dp 



ln/(z|p) 



dp 



ln/(p) 



(4) 



(5) 



(6) 



represent the contributions originating from the observation 
vector (i.e., the so called conditional information) and that due 
to a priori information, respectively. The last decomposition 
allows us to analyse the two above mentioned sources of 
information in a separate fashion. Therefore, in the following, 
we mainly focus on the evaluation of the matrix J p dSl, since 
this term unveils the contribution of map information in the 



estimation of agent position. As already mentioned in Sec. |II-A 
a smoothed uniform map model is adopted in the derivation of 
J p (analytical details are provided in Appendix [A]) so that its 
pdf /(p), unlike its (discontinuous) uniform counterpart ([TJ, 
satisfies the regularity condition p0| 



d_ 
dp 



in/(p; 







(7) 



(3) 



required for the evaluation of the Bayesian Fisher information. 

Moreover it is assumed that: 

1) the smoothing of /(p) is described by a smoothing 
function s(t), which is required to be a continuous and 
differentiable pdf for which the a priori Fisher informa- 
tion (FI) J s = E t I ( dy* ) \ can be evaluated; 
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2) the smoothing of /(p) along the x axis is independent 
from that along the y axis. 

In Appendix [A] it is proved that, given the assumptions 
illustrated above, J p d6l) can be approximated as 



diag 



E 



dy 



E 



dx 



(8) 

Note that such an expression is independent from the observa- 
tion model and depends from the smoothing function through 
its FI J s onlj|^] Given J p , the computation of the BCRB 
requires the knowledge of J z | p |5jl (see (|4jl); for its evaluation, 
instead of considering a specific observation model, we just 
assume that the cross-information about x and y is negligible, 
so that such a matrix can be put in the forrrj^] 



J Z | P — diag { Jjsja;, J Z |y } 



(9) 



Then, the BCRB associated with the pdf /(p) can be put in 
the form (see ([HJ, Q and 



e y b 




(10) 

It can be shown that the last two expressions are exact, and 
can be put in a closed form (not involving any integral), if 1Z 
is the union of a set of disjoint rectangles (all having parallel 
sides). 

From (|8]l and ( fT0] i it can be inferred that: 

1) Map-aware estimators should be expected to be always 
more accurate than their map-unaware counterparts; in 
fact J p dHJis positive definite, so that the trace of the 
matrix B (|10|> (which bounds the accuracy of map-aware 
estimators) is smaller than the trace of 



p)' 



diag 



{ J z|l' J z|y} 



(see (|9]l) which represents the Cramer-Rao bound (CRB) 
in the considered scenario (and bounds the accuracy of 
map-unaware estimators). Note that this is a well-known 
result (the contribution of a priori information is positive 
and lowers the BCRB), but may be concealed by the 



complex analytical form of (lOi 



2) 



The BCRB (lOi exhibits a complicated dependence on 
the map properties and, in particular, a complex non 
linear dependence on the map area An\ in fact, this 
parameter can be related to the functions {h m (x)} and 



In Section 



III-B 



the value to be assigned to J s is derived for the case of 
rectangular maps. 

3 The following expression becomes J z | p = diag | cr^ 1 , cry 1 1 if the 
specific observation model J2j is adopted; this result will be exploited later 
for comparison with other bounds derived for the model J2J. 



{ w n(y)} appearing in the right-hand side (RHS) of ( lOi 

as 



An= J 2^ w n {y)dy 



X 



E 



h m (x)dx 



However, the following asymptotic result suggests that 
the BCRB is expected to increase (and thus the accuracy 
of map-aware estimation to worsen) as An gets larger: 
if An — > oo the BCRB tends to the CRB (since the map 
support 1Z — > M 2 and the contribution coming from prior 
information vanishes^]) which always takes on a bigger 
value (see point 1). 

3) Only the smallest elements in the set {w n (y)} 
({h m (x)}) significantly contribute to the sum 
^ne^(y) ^(y) (EmeA^z) T^ij) appearing in (g) 
and ( fTO] ) and, consequently, may appreciably influence 
estimation accuracy; in the limit, if w n (y) — > 
(h m (x) — > 0) for some n and y (to and x), then the 
information about x (y) coordinate becomes infinite 
because the agent is constrained to lie, for those values 
of n and y (to and x), in a single point of the map 
support. 

4) The matrices J p ([8]l and B dTOj can be easily put in 
a closed form when {ft, m (x)}and {w n (y)} are step 
functions (i.e., the map is made by union of rectangles). 

Further insights can be obtained taking a specific class of maps 
into consideration. In particular, here we focus on the case of 
a map without obstructions, i.e., characterized by a support 7Z 
such that Afh(y) = and Afy(x) — \/x, y (note that convex 
maps belong to this class); under these assumptions, it is easy 
to show that J p dSl can be simplified as 



A 



n 



diag 



dy 



y wi(y)jx h i( x ) 



dx 



(11) 



From the RHS of the last expression it can be easily inferred 
that: 

1) Even in this case the dependence of J p from A-r. is 
complicated. However, analysing ( flT) , it is easy to note 
that if 1Z — > R 2 (so that An — > oo) the elements of 
J p tend to zero, so that the impact of prior information 
vanishes, as already mentioned above. 

2) The amount of information concerning the agent po- 
sition along the x (y) axis direction is proportional 
to J y w(y)dy (J x h(x)dx), where w(y) = l/w x {y) 
(h(x) = l/h\{x)); the last integral can be interpreted 
as the area of a "virtual map" characterized by the same 
projection y (X) as the real map, but whose width 
(height) at each point yE y (xE X) is the reciprocal 
of the width (height) of the real map. 

3) When prior information dominates over conditional in- 
formation, the map support 1Z minimizing the trace of 
the BCRB under the constraint of a constant area An 
has a square shape. In fact, if 1Z is a rectangle whose 
sides have lengths L x and L y , it is easy to show that 



4 if n - 

4.2]). 



the prior expressed by Q becomes "improper" (see |2I | Sec. 
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E >: Z = diag 



2A 



n Jy 



2A K 



x 



E c 
E c 

mGjV°(a;) 



w n (y) 



hm (^) 



y E Co 
E c 



Aw„(y) w n -i(y) w n+1 (y) 



Ah m (x) h m -i(x) h m+ i(x) 

5 ) 

dii (Tii (Jii 



dy, 



dx 



(14) 



we set L T 



k dia § { t ■ } = J * dia e { 4 > xr I ; then ' if 

7-L, L y = ^, where 7 and i denote two 
real positive parameters, the trace of the inverse of ( [TT] ) 
can be put in the form k- {^) 2 + I/7 2 ) and is minimized 
for 7 = 1 (i.e., for L x = L y ). 
4) In (jTTJ the only contribution coming from the smoothing 
function s(t) is represented by the factor J s . 
Further comments about the BCRB ( [T0| will be provided in 
Section IV where this bound is computed for some specific 
maps. 

Finally, it is important to point out that: a) in localization 
problems BCRB's often provide useful insights J9), fTO), fT2), 
but these bounds are usually loose for low SNR conditions 
[22 1, [23], i.e., when a priori information (the map in this 
context) plays a critical role due to the poor quality of 
observations; b) the BCRB analysis requires the adoption of 
the smoothed uniform pdf model for prior information (see 
Section |II-A| i. 

These considerations motivate the search for other bounds 
and, in particular, for the EZZB, which is usually tighter than 
the BCRB at low SNRs |23|- (27| a nd does not require Q to 
hold (so that the uniform model (fill can be employed as it is); 
this is further discussed in the following section. 

B. Zik-Zakai Bounds 

Similarly to the EZZB for 2D localization can be 
expressed as 

E 4 diag {E x , E y } >"L = diag {Z x , Z y } (12) 
where (see Appendix |Bj) 

Z --\JJ [f(P) + f(P + he -)\ p l (P, P + he,) h dpdh 

(13) 

v € {x,y}, e x 4 [1,0] T , e y 4 [0, 1] T , P^p,p + he u ) 
represents the error probability referring to the likelihood ratio 
test in a binary detection problem involving the hypotheses 

H : P = P 

H 1 : p = p + he v 

and the observation z; finally, the integration domain is 

V v = {(h, p) ■ h > A / (p) > A / (p + he v ) > 0} C M 3 

The EZZB (fT2]>-( 13 1 deserves the following comments: 
1) its formulation is obtained from the standard formulation 



(see 1 25 Eq. (32)]) removing the so-called valley -filling 
function and the maximization operator, as shown in 



Appendix [B] these modifications simplify the derivation 
of the bound at the price of a small reduction in its 
tightness; 

2) it exhibits a complicated dependence on the observation 
model, which comes into play in the evaluation of the 
probability P^ n (p, p + he u )\ 

3) it cannot be put in a form similar to Q, so that, generally 
speaking, the contribution coming from conditional in- 
formation (i.e., observations) cannot be easily separated 
from that associated with a priori information (map 
awareness in this case). 

As shown in Appendix [B] for a uniform map and the obser- 
vation model (|2|l the EZZB ( 13 1 can be put in the form ( 14 1 
shown at the top of this page where 



C(p) — / (p~ u ) u er f c 



( — 



du 



(15) 



Cov (fiA,Pl,P2) — 

rPA+P2 

(u-p A )uerfc . 

'PA \2VZ 



du- 



P2 u erfc 



PA+P2 



2V2 



du- 



PA+P1+P2 

(PA + Pi + P2 - w)uerfc . 
pa+pi \2v2 



du (16) 



An accurate analysis of the EZZB 
the BCRB 



14 1 and a comparison with 



10 1 lead to the following conclusions: 

1) The EZZB (fT4] ) exhibits a similar dependence as 
the BCRB (|10| in terms of the quantities {w n (y)}, 
{h m (y)}. In fact, if {w n (y)} and/or {h m (x)} increases, 
the trace of matrix Z in ( fT4| gets larger, since ( ( 15 1 is 
a monotonically increasing function. 

2) In the same environment the EZZB ( p"4} takes on larger 
values than the BCRB ( fT0) i and results to be a bound 
tighter to the performance of optimal estimators. This is 
due to the fact that the terms ( appearing in the RHS of 
( fl"4| ) give the EZZB approximately the same behaviour 
of the BCRB; however, the EZZB also contains terms 
based on C ou ( [To} , which is a positive valued function. 

3) The EZZB (fT^unlike the BCRB, is influenced by the 
shape of map obstructions (through the ( ov function and 
the {Aw n (y)}, {Ah m (y)} parameters). 

4) The two terms of ( fl4| depending on the function £ ov 
( 16 1 are influenced by the spatial extensions {Aw n (y)} 
and {Ah m (y)} of the obstructions separating consec- 
utive map segments (which are characterized by the 
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lengths {w n -i(y),w n+ i(y)} and {/i m _i(y), h m +i{y)}, 
respectively); in particular, the impact of obstructions is 
significant when Aw n (y) (Ah m (y)) is greater than one 
of {w n ^ 1 (y),w n+1 (y)} ({h m -i(y),h m+1 (y)}). 
5) The integrals in the RHS of ( p"4| ) can be easily put in 
a closed form when {h m (x)} and {w n (y)} are step 
functions (i.e., the map is made by union of rectangles). 
Similarly to the BCRB, further insights can be obtained 
considering the class of maps without obstructions. In this 



case the EZZB (14i simplifies as 



E >- Z 



2A 



diag 



c 



,V 



wi{y) 



hi(x) 



dy, 



dx 



(17) 



It is also worth mentioning that, if the map support 1Z can 
be expressed as the union of a set of disjoint rectangles (all 
having parallel sides), the EZZB ( [17} can be simplified further, 
obtaining a closed-form bound. 



The result ( 17 1 deserves the following comments: 
1) The structure of its RHS is similar to that of the 

Despite the differences between 



equivalent BFIM ( 1 1 



the integrand functions, our numerical results show that 
in this case both the EZZB and the BCRB are tight 
and that the EZZB does not provide additional hints 
on localization accuracy with respect to the BCRB (see 
Section |IV} , 

2) The EZZB predicts the correct variance of map-aware 
estimation when the SNR is very low whereas a specific 
assumption about the value of J s is required to obtain 
the same result in the BCRB case. For instance, let us 
consider a uniform map whose support 1Z is a rectangle 
whose sides have lengths L x and L y ; with model (|2]l and 
for <j x — } oo and a y — > oo (i.e., for a SNR approaching 
zero), the BMSE of the minimum mean square error 
(MMSE) estimator tends to diag {L 2 X /12, L 2 y /12). The 
EZZB provides this exact result since the RHS of ( [T7j ) 
can be put in the closed form 



Z 



and ((p)/(2p) 




(18) 



/12 for a SNR p approaching 



zerc[^] This result can be compared with its BCRB 
counterpart; in this case, from ( fTI) it can be inferred 
that diag{L|/ J s , L%/ J s }, so that J s = 12 should 
be selected to match the correct variance. 



3) The EZZB (17i confirms that map-aware estimators 



should be expected to be more accurate than their map- 
unaware counterparts, at least for maps without obstruc- 
tions: for instance, if a rectangular map of sides L x , L y 
is considered, then (TTSj is easily proved to be always 
smaller than diag { o^a y } (because < Q{p)/(2p) < 
1) which is the trivial EZZB^jfor the observation model 

5 This proof requires developing a Taylor expansion of the function 
<(p)/(2p) around p = 0. 

6 This result is obtained from (liil after proving that £(p)/(2p) — ■> 1 for 
p -> +oo. 



Q when the map is unavailable (i.e., when no prior 
information is available). This derives from the fact that 
a priori information lowers the EZZB (similarly to what 
occurs in the case of the CRB and the BCRB); this is a 



well-known result, but may not be evident in ( 17 I. 
4) Similarly to the BCRB, the dependence on the map area 
is complicated; however, it should be expected that the 
elements of Z get larger when An increases as sugges- 
ted by the following asymptotic result: when A-r — > oo 
the contribution coming from a priori information van- 
ishes and the EZZB dlSb tends to diag {a x , tr^}, i.e., to 
the map-unaware EZZB, which is larger (see point 3). 
As already mentioned above, the EZZB allows us to under- 
stand the role played by some map features, which do not 
have any impact on the BCRB. However, our numerical results 
evidence that, for complex maps and, in particular, at very 
low SNRs, this bound turns out to be somewhat far from the 
performance offered by a MMSE estimator. This has motivated 
the derivation of the WWB |28) , which can be tighter than the 
EZZB at very low-SNRs. 



C. Wess-Weinstein Bounds 

Similarly to the BCRB and the EZZB (see Q and ( [T2| , 
respectively), the WWB, for 2-D localization, can be expressed 
as 

E >r W = diag {W x , W y } 

where (29) 

(hA,p{L* (z;p+,p)} W 



(19) 



W v = sup 



z.p 



Ll (z;p+,p) - L\ (z;p , p) 



(20) 



and 



+ ^ 
P = P 



p = p - h v e„ 



are functions of the h u parameter; moreover 

A /(z,Pl) /(z|Pl)/(Pl) 



L(z;pi,p 2 ) 



(21) 



/(Z,P 2 ) /(z|p 2 )/(P2) 

is a likelihood ratio and v € {x, y}. Note that ( p0) > can be ob- 
tained from J29] Eqs. (7)-(8)] setting N = 2, s x = s 2 = \ for 
the WWB optimization parameter^] and H = diag {h x , h y }, 
since we are interested in a bound on the variance (not on the 
co-variance) of the estimates of p. 

As shown in Appendix [C] for a uniform map and the 
observation model (|2]i the WWB ( 20 1 can be put in the form 
(22 1 shown at the top of the next page, where h x , h y denote 



the parameters to be optimized to make the bound tighter, 

X x (h x ,TZ)= / I k (p)Ik(p + h x e x )dp 



= / ^2 uj(w n (y),h x )dy+ 

I ^2 UJ ov{^w n (y),w n ^ 1 (y),w n+ i(y),h x )dy (23) 

7 These choices make the bound tighter, as pointed out in [28] . 
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E y W = diag 



sup 







X x [hxi 




2 


2A K 


X x (h x ,Tl) - exp 


(-m 


lx (h x ,TZ) 



, sup 

h v £R 





h v ex p(-^|) 


Xy (hy, 


n) 


i 


2An 


X y (h y ,lZ) — exp 


(-3) 





(22) 



and 



= 2 f V Lj(lw n (y),h x )dy (24) 

are functions of the parameters h x , h y and of the geometrical 
features of the map support TZ (dual expressions hold for 
X y (h y ,lZ) and j y (h y ,TZ)); finally 



uj {w, h x ) = (w - h x )l[ 0;w ] (h x ) 



(25) 



oj ov (Aw,w 1 ,w 2 ,h) = 

h — Aw h G [Aw; Aw + w 2 ] 
w 2 h e [Aw + w 2 ; Aw + u>i] (26) 

w tot -h h G [Aw + wi\ Wtot] 

w tot = Aw + w x +w 2 ■ 



Note that the structure of (h x ,1Z) ( 23 i is similar to that of 
some terms contained in the RHS of the EZZB ( [14) ; however, 
in the WWB the noise variances appear in some exponential 
functions (f22]i only. 



Unfortunately, it is hard to infer from ( 22 1 any conclusion 
about the dependence of the WWB on the features of a given 
map; this is mainly due to the presence of the sup (■) operators 
and to the discontinuous behaviour of the functions to (•) and 
uj ov (■). Despite this, the WWB may play an important role, 



since, as shown in Section IV for specific maps, it is tighter 
than the BCRB and the EZZB at low-SNRs. 

In the following Section the general bounds provided above 
are evaluated for specific maps, so that some additional 
insights about the ultimate accuracy achieved by map-aware 
estimators can be obtained. 

IV. Numerical Results 

The bounds illustrated in the previous Section have been 
evaluated for the following 3 specific maps: 

1) map #1 : a unidimensional map whose support 1Z consists 
of N r = 2 disjoint segments (denoted IZi and 1Z 2 in the 
following) having the same length w and separated by 
Ax meters; 

2) map #2: the "L-shaped" 2-D map shown in Fig. [2] 
and fully described by the set {wi,w 2 , hi,h 2 } of 
geometrical parameters; 

3) map #3: the 2-D map shown in Fig. [3] note that such 
a map models the floor of a large building and is quite 
complex, so that its geometrical parameters are not given 
and analytical bounds are not evaluated for it; instead, 





w 2 


* * 















X 



Figure 2. Support region 7Z of a "L-shaped" map (map #2). 




X 



Figure 3. Support region TZ of a map modelling the floor of a big building 
(map #3). 



some numerical results generated by means of computer 
simulations will be illustrated. 

It is important to point out that maps #2 and #3 are ex- 
amples of 2-D rectangular uniform maps, i.e., each of them 
is characterised by a support TZ that can be represented as 
the union of non-overlapping rectangles, all having parallel 
sides; such non-overlapping rectangles form the so called 
rectangular covering of the map support TZ |30|-|32). Note 
that any map characterised by a support of arbitrary shape 
can be always approximated, up to an arbitrary precision, by a 
rectangular map if a proper covering, containing a sufficiently 
large number of rectangles, is selected. For this reason, maps 
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#2 and #3 are relevant and practical examples. Map #1 can be 
also deemed relevant from a practical point of view since it can 
be thought of as a transversal slice of a 2-D map (consisting, 
for instance, of two road lanes or two rooms facing each other); 
in addition, such a simple scenario allows to acquire more 
insights than maps #2 and #3, as shown later. 

In all cases the observation model ([2]) has been adopted 
(with X = a 2 and X = <r 2 I 2 in the 1-D and 2-D scenarios, 
respectively) and the tightness of the bounds has been assessed 
comparing them with the root mean square error (RMSE) 
performance of optimal estimators evaluated by means of 
extensive computer simulations^] In addition, in accordance 
with the theoretical arguments of Paragraph III-B (see the 



w = 1 m, Ax = 1 m 



comments to (fT7]i), J s = 12 has been selected in the evaluation 
of all the BCRB results illustrated in this Section. Finally, 
in our simulations N runs = 10 4 realizations of the random 
parameter p have been generated according to a uniform 
distribution over the map support TZ; then Gaussian noise 
samples (characterized by a — 3 m if not explicitly stated) 
have been added to generate the noisy observations, according 
to 0. 

A. Map #1 

Bounds - The BCRB, the EZZB and WWB for map #1 are 



(see (10 1, (14 1, and (22 1, respectively) 



E x > B x = o l 



4p 



[2C(P) + C™(PA,P,P)] 



(27) 
(28) 



and E x > W x , respectively, where 

h 2 exp(-^) A 2 (h,w,Ax) 



W x = sup 



hem. Aw [A (h, w, Ax) — exp (—5^2) j(h,w,A 



(29) 



p = w/a is the SNR, p& = Ax /a, 

A (h, w, Ax) A 2uj (w, h) + uj ov (Ax, w, w, h) 

and j(h,w,Ax) A 2uj(hw, h). It is important to point 
out that that no approximations have been adopted in the 
derivations of the BCRB, the EZZB and the WWB for the 
considered 1-D scenario, so that (|27|-([29| represent exact 
bounds (of course, a smoothed uniform map has been assumed 
for the BCRB only). 

Estimators - Map-aware localization can be accomplished 
exploiting the MMSE estimator 



^mmse(z) = (l/c(z)) V] / xAf(z;x, a 2 ) dx 



(30) 



where c(z) A J2 n =i Itz N ( z ! x > (j2 ) dx is a normalization 
factor, or the maximum-a-posteriori (MAP) estimator 

r 2\ 



Xmap{z) — argmaxjV (z; x, a 2 ) 



(31) 



8 The software developed to generate our numerical results is available at 
http://frm.users.sf.net/publications.html in accordance with the 
philosophy of reproducible research standard 1331. 
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Figure 4. RMSE versus the SNR p for the MMSE {30}, MAP {31} and ML 
estimators. Map #1, w = 1 m and Ax = 1 m are assumed. The BCRB \27\ , 
EZZB {28} and WWB {29} are shown for comparison. 
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Figure 5. RMSE versus Ax for the MMSE {30}, MAP (IT} and ML 
estimators. Map #1, w = la and a = 3m are assumed. The BCRB (27} , 
EZZB {28} and WWB {29} are shown for comparison. 



On the contrary, if the map is unknown, the (map-unaware) 
maximum likelihood (ML) estimator x ML (z) = z can be used 
(obviously, the resulting estimate has the same variance a 2 as 
the observation noise). Note that the RMSE performance of 
these estimators is influenced by the geometrical parameters 
w and Ax and by the noise standard deviation a. 
Numerical results - Figs. |4] and [5] compare the square roots 
of the BCRB, EZZB and WWB bounds y/B^, v%^, VWx 
and the RMSE of ^mmse^), x map (z) and x ML (z) versus p (with 
w and Aa; fixed) and Ax (with w and a fixed), respectively. 
These results evidence that: 

1) The BCRB is tight at high SNRs only; in such conditions 
all the bounds and the RMSE of both map-aware and 
map-unaware estimators tend to a (see Fig. |4|: condi- 
tional knowledge dominates. 
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RMSE for the MMSE estimator on map #1 




Figure 6. RMSE versus Ax and w for the MMSE estimator pO) . Map #1 
and a = 3 m are assumed. The data has been smoothed to remove simulation 
noise and improve readability. 



2) Unlike the BCRB, the EZZB and, in particular, the 
WWB exhibit a dependence on Ax and are much tighter 
(see Fig. |5]l to the performance of optimal estimators 
and, for this reason, represent more useful tools for 
predicting system accuracy. 

3) The role played by map knowledge in estimation ac- 
curacy becomes significant as p decreases, i.e., as the 
quality of the observations gets worse; indeed map 
knowledge prevents the bounds and the RMSE of the 
map-aware estimators from diverging (note that in Fig. 
|4]the map-unaware RMSE diverges for p — > 0). 

4) The performance gap among the MMSE, the MAP and 
especially the ML estimators is significant for low SNRs 
(see Fig. [4}. 

5) If the gap Ax is small (i.e., if the segments IZi and 
IZ2 are close), the RMSE of the considered estimators 
is low. On the contrary, if Ai increases, the RMSE gets 
larger, reaches a maximum and then decreases, tending 
to w/\/l2 (see Fig. £5j. This behaviour can be explained 
as follows: if the noise level is large with respect to 
Ax, the estimate of the agent position may belong 
to the wrong segment, thus resulting in a large error. 
However, a further increase of Ax entails a reduction 
of the RMSE, since, for a large spacing between IZi and 
7Z2, it is unlikely that the wrong segment is selected in 
map-aware estimation. 

Further insights are provided by Fig. [6] which shows the 
RMSE versus w and Ax for the MMSE estimator (note that 
the MMSE data shown in Fig. [5] can be extracted from Fig. [6] 
setting w — 1 m). In fact, the numerical results shown in this 
figure evidence that: 

1) For w > 3cr the RMSE decreases only mildly as Ax 
gets larger and the main geometric feature of the map 
affecting estimation accuracy is the area of 1Z; in this 
case a map-aware estimator still provides an advantage 
over map-unaware counterparts (the maximum RMSE 
for the first estimator is 2.6 m, as shown in Fig. [6] and 



is smaller than a = 3 m, which represents the RMSE 
for the ML estimator, as shown in Fig. |5J. 
2) For w < 3fi the RMSE reaches its maximum for 
Ax ~ 3er — w; for this reason, any combination of 
the geometric parameters and noise level satisfying such 
equality should be avoided. 
Based on the numerical results shown in Fig. [6] we can 
draw some observations about an indoor localization system 
operating in two adjacent rooms, each having width w room 
and separated by a wall having thickness Ax wa u — 0.2 m, 
and processing noisy observations characterized by a = 3 m 
(realistic value for RSS systems, see [34|-[36|). From Fig. [6] 
we can infer that, even if w room < 3d = 9m, since Ax wa u 
is appreciably smaller than 3cr — w, the advantage in terms 
of RMSE provided by a map modelling wall obstructions 
is small with respect to that offered by a map modelling 
the environment as a single connected room having width 
%w room + Ax wa ii- In fact, in Fig. [6] if the presence of the 
wall is accounted for, the map-aware estimator performance 
is associated with the point (Ax wa iu w room ), whereas, if it is 
not, with the point (Q,w room + ^Ax wa u). Unless the room 
size w room is very small (i.e., w room — Ax wa u) the two 
"operating points" will be very close (due to the small value 
of Ax wa u) and the resulting RMSE values will be similar so 
that, in this case, a map modelling wall obstructions is not very 
useful. Note also that, if the quality of observations improves 
substantially (i.e., a — > 0), map information plays a less signi- 
ficant role than the conditional information provided by such 
observations, so that a detailed a priori map modelling walls is 
not useful in this case, too. On the contrary, maps may play an 
important role in localization in outdoor environments, where 
the region characterized by w < 3cr and Ax <C 3er — w in the 
(Ax, w) plane is easier to reach, since noisy observations and 
large obstructions should be expected; for this reason, in this 
case the presence of the obstructions in map modelling should 
not be neglected. 

B. Map #2 

Bounds - Let us now focus on map #2. The BCRB and the 
EZZB in this case are (see (JTOj and ( [T4| , respectively) 



E h B = diag { a 1 



p\ 



Pa 



( hx 

\V>1 



W-L+W2 



Pa 



" 2 - J, 



tl'2 

hi 



E " Z = diag {^ 

a 2 
2PA 



PA "T \ /, ; T hx+h 
W\ + W2 ' 



(32) 



WiC 



17 

hi + h 
a 



h 2 ( 



+ w 2 ( 



a 
hi 
a 



(33) 



respectively, where p_^ = A-jz/a. The WWB derived from 



22|) is shown in (|34|) at the top of the next page; in ( 34 1 the 



\ x (h x ,TZ) and j x (h x ,lZ) functions have the expressions 
\ x (h Xl TZ) = hiu(wi + W2, h x ) + h2U>(wi,h x ) 



10 
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E y W = diag < sup 





h 2 x exp ( 


K; fox. 


n) 


2An 


X x (h x ,K) 




j x (h x ,TZ) 



sup 





h 2 y exp ( 


ia 2 ) 


xl (h y ,n) 


2An 


\ y (h y ,K) 


— exp 


( 



(34) 



7z (h x ,K) = 2/iiw 



/ Wi + w 2 



2h 2 u ( 



Wl 

2 ' 



wi = 5 m, h-2 = 5 m, tr = 3 m 



respectively (dual expressions hold for \ y (h y ,lZ) and 
7 y (h y ,lZ)). Note that hi and /12 (wi and u^) play a symmetric 
role in (32) , p3j ) and (34i; this suggests that a change in the 
geometrical features improving B x will reduce B y and vice 
versa, when the overall area An of the map is kept constant. 
Estimators - Map-aware localization can be accomplished 
exploiting the MMSE estimator 



e(z) 



*^MMSE 

2/mmse 



= (l/c(z)) // pAA(z;p,S) 



n 



dp 

(35) 

where c(z) = JJ" K A/" (z; p, S) dp. On the contrary, the (map- 
unaware) ML estimate is given by p M l(z) = z- 
Numerical results - Fig. [7] shows the square roots of the 
BCRB components \JB X , ^f~B~ y , its trace y / tr{B} and the 
RMSE of the estimators a; MMSE , 2/ M mse> Pmmse 

(z) and p ML (z) (the 



EZZB (33 1 and WWB (34i are not shown to ease the reading); 
a = 3m, wi = 5m and h 2 = 5m have been assumed in this 
case. In addition, when increasing w 2 , hi is reduced according 
to the law 

, / x An - h 2 wi 

hi{w 2 ) = ■ , 

Wi + w 2 

so that the equality An = 75 m 2 always holds. From these 
results it is inferred that: 

1) As already noted for map #1, map-unaware estimation, 
which is unaffected by a change in the geometrical 
features of the map, provides a worse accuracy than its 
map-aware counterpart (the RMSE of ML estimator is 
~ 21% larger than that of its MMSE counterpart). 

2) Increasing w 2 and decreasing hi changes the aspect ratio 
of a portion of the map (see Fig. [2]); this results in an 
appreciable reduction in the RMSE of the parameter y, 
but also in a mild increase of the RMSE of the parameter 
x. Therefore, the improvement along one direction is 
compensated for by the worsening in the orthogonal 
direction, so that no significant advantage for the RMSE 
of 

Pmmse ( z ) IS found. 

The last result suggests that when the accuracy in the estim- 
ation of a given element of p is improved maintaining An 
constant, a reduction in the estimation accuracy of the other 
element should be expected, so that that the overall accuracy 
does not change much. 

C. Map #3 

Observation model - Finally, let us analyse the RMSE per- 
formance of localization systems operating in the environment 
described by map #3 in two cases: A) the map-aware estimator 
knows map #3 exactly; B) the map-aware estimator assumes 
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Figure 7. RMSE versus w 2 for the MMSE estimator |35j. Map #2 (Fig. 
|5J, a = 3 m and A-jz = 75 m 2 are assumed. The BCRB components y/B x , 
y/By, ^/tr{B} (see j32J) and the RMSE for the ML estimator are also 
shown for comparison. 



that the localization environment is described by the bounding 
box of map #3. Note that the last case is significant from 
a practical point of view, since, in the absence of detailed 
information about the floor of a building, the bounding box 
represents the best a priori model which can be adopted to 
describe the floor itself. 

In addition, in this Paragraph model Q is compared against a 
more specific counterpart, called "ranging observation model" 
in the following; for such modelling we assume that 1) the 
observations are acquired by 4 anchors placed in the map 
corners; 2) the observation provided by the i-th anchor can 
be expressed a^] 



Zi = Hp -Pi 



(36) 



with i = 1,2,3,4, where rij represents Gaussian noise with 
variance a 1 and pf represents the i-th anchor position. 

The scenario considered in this example is helpful to assess 
1) to what extent a detailed knowledge of the geometrical 
features of a map can improve localization accuracy with 
respect to the availability of a simplified map model (like a 
bounding box model) and 2) if a localization system based 
on the Gaussian model (|2]) can exhibit a similar behaviour as 
that of a localization system based on ( 36 1, at least in the 
considered cases. 

Estimators - Given the model (f36l), the MAP estimator can 



This model is commonly used in the study of localization systems 
exploiting TOA and, with some slight changes, RSS. 
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be expressed as 



fMAP 



(z) 



axgmaxjVfz;/i(p),S fl> ) (37) 



where M (p) 4 [||p _ p?|| , ||p - p|||] J and S H = a 2 I 4 . 
If the observation model Q is adopted in place of ( [36| , the 
MAP estimator 



fMAP 



(z) = arg max 

-,G 



Af(z;p, 



(38) 



is obtained, where X = S = a 2 l2- 
Numerical results - Fig. [8] compares the RMSE performance 
achieved by the MAP estimators (37i and (38 i versus a in 
cases A and B. These results show that: 

1) Surprisingly the selection of the map model (map #3 in 
case A or its bounding box in case B) does not affect the 
MMSE estimation performance, whatever the model for 
the observations. This result can be motivated as follows. 
For each value of a, there is some portion of the map 
where the combination of noise level and geometrical 
features (see the analysis of map #1) entails a large 
RMSE and where a detailed map (map #3) does not 
provide a substantial advantage over simpler maps (e.g., 
the bounding box); since the Bayesian RMSE results 
from an average of the agent position over the entire 
map, it is biased by such portions of the map for each 
value of a. Such an effect influences also any global 
Bayesian bound, like the BCRB, the EZZB and the 
WWB: their values result from an average over the entire 
parameter space (in this case 1Z) and can be strongly 
influenced by the parameter values (in this case, the 
values of p) characterized by the largest errors [37) . 
Note also that local bounds (e.g., CRB), which model 
the parameters to be estimated as deterministic variables, 
do not suffer from this problem, but cannot exploit 
a priori information and represent performance limits 
for unbiased estimators only [25]. In our case study a 
detailed knowledge of map #3 (case A) provides, in 
various subsets of 1Z, a substantial improvement with 
respect to the knowledge of the bounding box (case B), 
but the effects of such subsets are then "averaged out" 
when evaluating the Bayesian RMSE, as shown in Fig. 
i 

2) The adoption of the Gaussian observation model |2]) 
entails a worse accuracy than that provided by the 



ranging observation model (36i. This can be related to 
the placement of the anchors and to their number; in 
fact, range observations coming from anchors displaced 
with a constant angular spacing allow to average out the 
measurement errors better than the positional observa- 
tion of the Gaussian model |2} when the same noise 
variance a 2 is assumed in both model£3 If the offset 

10 It is important to note that using the same number of observations for both 
models may generate misleading results. In fact, on the one hand the minimal 
number of measurements usually considered for the ranging observation model 
is 3 (this allows non-ambiguous 2-D localization); our choice of 4 anchors 
and 4 measurements is perhaps more realistic. On the other hand, for the 
Gaussian observation model, 2 positional observations (one about x and one 
about y) are considered since 2-D localization is performed. 
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Figure 8. RMSE versus a for the MAP estimators j37| and (38) , using either 
the Gaussian observation model (see J5J) or the ranging observation model 
(see (36) ). The a priori map is map #3 (case A) or its bounding box (case 
B). The RMSE of ML estimators is also shown for comparison. 



between the curves referring to the Gaussian and the 
ranging model is neglected, our results show that both 
curves exhibit the same behaviour, when a is varied; this 
is an hint to the fact that the insights obtained from the 
bounds derived for the model Q may be useful also in 
the analysis of systems employing different observation 
models (like (f36|). 
Finally, comment 1) suggests that, when assessing the im- 
provement in localization accuracy coming from map aware- 
ness, the evaluation of the overall RMSE does not provide 
sufficient information about the importance of a given map. 
In fact, if the effects of different portions of the maps are 
analysed, it is found that some of them play a more relevant 
role than other ones; in this perspective, the insights and the 
intuition acquired from the results referring to map #1 are very 
useful, given that map #1 can be thought of as a transversal 
slice of a subset of any 2-D map. 

V. Conclusion 

In this manuscript the impact of map awareness on lo- 
calization performance has been investigated from a theoret- 
ical perspective, evaluating different accuracy bounds. Such 
bounds provide some general indications about the role and 
importance of this form of a priori information. Our study has 
evidenced that, unluckily, the tighter is an accuracy bound, 
the more complicated is its dependence on the geometric 
parameters of maps. This has motivated the analysis of the 
developed bounds in specific environments; in particular, three 
different maps have been considered and, for each of them, 
the accuracy bounds have been evaluated and compared with 
the performance offered by map-aware and map-unaware es- 
timators. Our analytical and numerical results have evidenced 
that: a) map-aware estimation accuracy can be related to some 
features of the map (e.g., its shape and area) even though, in 
general, the relation is complicated; b) maps are really useful 
in the presence of some combination of low SNRs and specific 
geometrical features of the map (e.g., the size of obstructions); 
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c) for some combinations of SNRs and geometrical features, 
there is no need of refined maps since knowledge of map 
details provides only a negligible performance gain; d) in a 
given environment the EZZB and, in particular, the WWB are 
usually much tighter than the BCRB; e) in the cases of maps 
containing many features, the importance of map awareness is 
better captured by the analysis of subsets of those map. 

Future improvements include the development of map- 
aware bounds for specific observation models; in particular 
models including bias due to the non-line of sight (NLOS) 
propagation may unveil important insights about map-aware 
localization systems operating in NLOS conditions. 

Finally, it is important to point out that all the bounds 
derived in this manuscript can be also exploited for other 
estimation problems (provided that the vector of random 
parameters to be estimated is characterized by a uniform 
a priori pdf) and can be easily generalized to the case of 
estimation of an TV -dimensional parameter vector, with N > 2. 

Appendix A 
Derivation of the BFIM and BCRB 



In this Appendix the derivation of d8b and ( 1 1 1 is sketched 



To begin, we consider a real function s(t) having the following 
properties: 

P.l) it is continuous and differentiable; 
P.2) s(i)>0VteM; 
P.3) / s(t)dt = 1; 
P.4) f^dt = 0; 

P.5) its support is the intervaj^] [-1/2; 1/2]; 
P.6) s(0) = 1. 

These properties have the following implications: a) s(t) is a 
pdf function (see P. 1-P.3); b) an a-priori FI J s can be computed 
for it (the regularity condition P.4 ensures the FI existence) 
[20 1; c) s(-) can be used to model bounded statistical distribu- 
tions (i.e., maps), eventually after some scaling and translation 
(see P.5, P.6). As far as the last point is concerned, it is worth 
mentioning that 

~s(t,a,b)±(l/b)-s((t-a)/b) , 

with b > 0, is still a pdf function and its FI is J s (b) — J s /b 2 . 
Any function s(-) sharing the properties P.l -P.6 represents a 
"smoothing function"; specific examples of smoothing func- 
tions can be found in J38] . 

Any smoothing function s(t) can be used to define the 
analytical model of a smoothed uniform map. To show how 
this can be done, we start considering a unidimensional (1-D) 
scenario first, where the agent position to estimate is the scalar 
x. The support 72. C M of a 1-D uniform map can be always 
represented as the union of N r disjoint segments spaced by 
N r — 1 segments where the map pdf f(x) is equal to 0. Let the 
N r segments of the support 1Z be indexed by the odd numbers 
of the set M° = {1, 3, 2N r — 1} and let c n (w n ) denote the 
centre (width) of the n-th segment, with n € Af°. Then the 

"The measurement unit for the support of s(t) is the same as that adopted 
for p and z (see Sec. |nj. 



map pdf f(x) associated with this scenario can be expressed 
as 

J_ ( x ~ Cn ^ 1 



n 



n€JV° 



W 



■R 



^ w n s(x,c n ,w n ) 



(39) 

where Wr = J2neAf° Wn - ^ ne a P r i° r i FI associated with 
f(x) is given by 



J T = E 




1 

1 



^2 w n J s (w n ) 




lns(x,c n ,w n ) 



n£JV"° 



E 



(40) 



Let us extend these results to 2-D maps. Exploiting some 
definitions given in Section|II]and assuming that the smoothing 
along x and that along y are independent, the pdf for an 
arbitrary smoothed uniform 2-D map can be put in the forrrp] 



Hp) 



i 



E - 
E - 

m£j\f°(x) 



w„(y) 
h m (x) 



(41) 



where c x>n (y) (c V}m (x)) and w n (y) (h m (x)) denote the centre 
and the length, respectively of the n-th (m-th) segment. Let 
us now compute the elements of the BFIM J p associated with 
this pdf (see (pi). The analysis of (41 1 can be connected to 



that of ( 39 1 exploiting the formula (iterated expectation) 




E, <^ E xly 



and adopting the approximation 

J_ \p ( x - c XiU (y) 



/(P) 



e n 

Wn(y) 1 
An Wr(v) 



(42) 



(43) 



E ■ 



w n {y) 



where Wn{y) - E„eAA°( y ) w n(y)- It is worth noting that in 
( 43 1 the smoothing along the y coordinate is ignored. This 
is a reasonable approximation whenever the smoothing only 
modifies a narrow area of the edges of 7Z compared to the area 
of the flat regions of ( |4T| ; the selection of a proper smoothing 
function ensures that such a condition holds (see Sec. |II-A| and 

mil). 



Thanks to (43 1, the evaluation of the inner expectation in (42 1 
is equivalent to that of the FI referring to a 1-D map which 
consists of N r = Nh{y) segments having widths {w n (y)} 

12 It can be easily proved that this pdf satisfies regularity condition 



E P {^)} 
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and centred around the points {c Xin (y)}; for this reason, from 
(|39l-(|40l it is easily inferred that 



E 



x\y 



'flln/(p) 
dx 



E 



i 



An Wn(y) ^ w n (y) 



(44) 



Substituting the last result in the RHS of (42i produces 

j s r i 



Mi i- ir I E —r^f(y) d y 



(45) 



where /(y) = J" f(p)dx is the pdf of j/ only and it has been 
ignored in the integral since the smoothing only affects a small 
portion of the integration domain (see Section [IIJ. A similar 
approach can be adopted in the evaluation of [J p ] 2 2 (ignoring, 
in this case, the smoothing along x); this leads to 



[J P ] 



x 



E 



hm (%) 



dx 



(46) 



Note also that the approximated results ( 45 1-( 46 1 become 
exact if TZ is a rectangle having its sides parallel to the 
reference axes, since in this case the smoothing along x is 
truly independent from the smoothing along y. For the same 



reason the cross terms 



[J p ] 2 x and 



'PJl,2 



are exactly equal to 



zero for a rectangle; however since 2-D maps usually exhibit 
a weak correlation between the variables x and y, these terms 
can be neglected; this leads to ([8]), from which (jTTJ is easily 
obtained assuming that Af%(y) = and Afy(x) — Vx,y. 

Appendix B 
Derivation of the EZZB 



In this Appendix the EZZB ( 13 1 is derived following a 
procedure conceptually similar to the one provided in [25] 
Sec. II.B]. We start from the identity (see B5] p. 24]) 



E v = ~j^Pr{\t u \>~\hdh (47) 

where = [p(z) — p] v . Following |25|, the probability 
appearing in ( |47] > can be put in the form 



Pr \\U> 



Pr <Cu > 



Pr 6< 



: J RN Pr > = Po } f( p o) dp o 
■jT^Pr|^<~|p = p 1 }/(p 1 )dp 1 



where /(•) denotes the pdf of p. Let p = p and p : = p + 
<5(/i), where is some function of the integration variable 
of ( |47j ). If the integration domain of ( |48| ) is constrained to be 
a subset 7Z(h) of the map support 1Z, such that /(p) > and 
f(p + S(h)) > 0, then, dividing and multiplying the RHS of 



<4jj) by /(p) + /(p 4 
eqs. (22)-(30)] 

Pr||6l>|l> 



■6(h)) > 0, produces the expression |25 



[/(p) + /(p + Pl(p, p + 5(/i))dp (49) 



(Pl(p,p + is defined in Sec. [IIFB| > which holds for 

any 6(h) satisfying the equality a 7 6(h) = h, where a £ M 2 . 
Note that ( |49] > generalizes p5] eq. (30)] to the case where the 
pdf /(■) of the parameter to estimate has a bounded support. 
Here a = e„ is selected so that 6(h) — he u . Then substituting 



(49 1 in (47 1 produces (13i, where the integration domains of 
h and p have been merged in the set V v . 



Like in the previous Appendix, let us evaluate now ( 13 1 for a 
1-D uniform map whose support KcR consists of the union 
of 7Y r disjoint segments. In the following such segments are in- 
dexed by the odd numbers of the set M° = {1, 3, 2N r - 1} 
and the lower (upper) limit of the segment n € N° is denoted 
l n — c n — -ff- ( u n — c n + 2 ^ L ), where c n (w n ) represents the 
centre (length) of the segment itself. Moreover, it is assumed, 
without any loss of generality, that c\ < C3 < ... < C2jv r -i- 
In this 1-D scenario ( fl3] > simplifies as 



1 



[/(r) + f{r + h)] P1(t, t + h)h dr dh (50) 



where p = r, 

V = {(t, h) : h > A /(r) > A /(r + fe) > 0} 

and P* in (r,T + h) represents the minimum error probability 
of a binary detector which, on the basis of a noisy datum z 
and a likelihood ratio test, has to select one of the following 
two hypotheses H : x = r and H\ : x — t + h. Note that, 
for any (r, h) € V, we have that /(r) + /(r + h) = 2/Wn 
(with Wr = X)nGAA° because of the uniformity of the 
considered map. If the Gaussian observation model |2]) adopted 
to the 1-D case is assumed, the average error probability of 
the above mentioned detector is easily shown to be 

p l(r, r + h) = ^ erfc (h/(2aV2] 

(where a denotes the standard deviation of the noise observa- 
tion model), which is independent of r; then, p0| ) simplifies 

as 



Z= J_ [[ hcTic (J^ 
2W K JJ V \2aV2 

= S., // uerfc ( — 
2Wn JJq \2y/2 



dr dh 
dt du 



(51) 



(48) where u = h/a, t = r/a and 



Q = {(t, u) : u > A f(ai) > A f(a(t + u)) > 0} 

can be shown to be 2-D domain consisting of the union of N r 
triangles and N r (N r — 1)/2 parallelograms in the plane (t, u), 
as exemplified by Fig. [9] which refers to the case N r — 3. In 
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particular, the contribution to the RHS of ( f5T] i from the i-th 
triangle 



Qi = 1 (t,u) :^<t<^A0<«<^-f 
a a a 



can be shown to be a 2 ((p.i)/(2pw), where pi = Wi/a, 
Pw — W-ft/cr, the function is defined in (15 1 and 



i G N° (see p8[ for more details). As far as the parallelograms 
are concerned, we retain only the contributions coming from 
(N r — 1) parallelograms 



(t,u) : - < t < — A- 
a a a 



t <u< 



with i £ Af e = {2,4, ...,2(N r - 1)} (see Fig. [9}; such con- 
tributions are given by a 2 ( ov (pA, Pi-i, p l+ i)/(2p w ), where 
p Ai = Axi/a, Axi = l t+ i - m-i and Cov(PA, Pi, P2) 
is defined in ( fT6) > for p\ > P2 (it can be shown that 
Cov (pa ,PuP2) = Cov (pa,P2,Pi))- Note that the contributions 
to pT| i coming from the discarded parallelograms are always 
positive since uerfc (y^p^j > for u > 0. Then, substituting 
the above mentioned contributions in ( 5 1 1 produces the lower 
bound 



Z > 



2p w 



)+ Cot; (PA n ,Pn-l)Pn+l) 



(52) 

for the EZZB. The approach developed for a 1-D uniform map 
can be extended to a 2-D scenario keeping into account that 
a) the term P^ la (p, P + he u ) of ( 13 1 takes on a similar form as 
P* m (r,T + h) of ( 50 1, b) the integral on p can be decomposed 



in a couple of nested integrals, one over X and the other one 
over y, and c) eq. ( |52] > can be exploited for the inner integral. 
Further details can be found in |38) . 

Appendix C 
Derivation of the WWB 

In this Appendix we first derive the WWB for a 2-D uniform 
map characterized by a support 1Z, assuming the observation 
model Q. Note that this derivation can be easily generalised 
to a A-D scenario. To begin, we note that, since /(z|p) = 
A/ r (z;p, S) (see (|2|i), the likelihood ratio _L(z;p!,p 2 ) (21 1 
is given by A/"(z; pi, S)/A/"(z; p 2 , S) Vpi, p 2 G 7vL Then, it 
is not difficult to prove that the expectation appearing in the 



numerator of ( 20 1 is given by 



E, 



.. p jiJ (z;p + h u e v ,p)\ 

^— [ [ A/"' (z;p + /i„e„,£)A/'5 (z;p,£)dpcfe 

^72. JR2 Jv v (h„) 



■ exp . 

-4k V oaf, j Jv„{K) 



dp 



(53) 



where V v (h v ) = {p : /(p) > A /(p + h u e v ) > 0} is a 
slice of the integration domain V v C M 3 involved in the 



evaluation of the EZZB [38] and v G {x,y}. The expectation 



appearing in the denominator of ( [19} can be expressed as 



E 



z.p 



L 2 (z; p + h v e v , p) - L 2 (z; p - /i„e„, p) 



= E Zi p {L (z; p + h u e v , p)} + E ZiP {£ (z; p - /i^e,,, p)} 
- 2 • E ZjP ji 2 " ( z; p + /i^e,,, p) (z; p - ft,„e„, p)| 

i.e., as the sum of three terms which are denoted A, B and 
C, respectively, in the following. It is important to note that 



A = E z , p {L (z; p + h u e v , p)} 

= -7— 1 [ Af(z; p + h v e v . H)dpdz 

■An 7k 2 Jp u (h u ) 

If 1 

= — / dp= —\ v {h vi n) 
■An Jv v (h v ) ■An 

= Ez.p {L (z; p - h„e u , p)} 

= B 



(54) 



and 



C = - 1 - f [ N{z;p + h„e v ,Y,)dpdz 

■An jr 2 Jv„(h„) 

hi 



1 

1 

An 



exp 



exp 



2^3 
2a?. 



dp 



(55) 



where 



V V {K) 4 {p : /(p) > A /(p + te^) > 
Nf{p-h,e v ) > 0} 

and the functions \ u (h v ,7V) and 7„ 7\L) are defined in 



(23 1 and d24J. Then, substituting d53Jl, d54b, (55 1 in (20 1 



produces (|22) 



Finally, it is worth pointing out that \ v (h v ,TV) (23i and 



^uihy^TV) (24 1 can be explicitly related to the functions 



{w n (-)}, {Aw n {-)}, {h m (-)}, {Ah rn (-)} introduced in Sec- 
tion [n] To clarify this point, let us focus on X x (h x ,lZ) and 
IxihxiTt) which can be derived, for a given ordinate y, 
analysing the integration domain 

V x (h x ,y) 4 {t : f(r, y) > A f(r + h x , y) > 0} . 

If the variables t = r/a and u = h x /a are defined, this 
domain can be put in the form 

Q x {u, y)={t: f{at, y) > A f(a(t + «), y) > 0} , 

which describes a slice of the (t, u) plane of Fig. 
[9] The contribution of the i-th triangle Qi to 

K(h x ,K(y)) ('fx{h x ,'R.{y))) is given by ui(wi(y),h x ) 
(u ov (Awi(y),Wi-i(y),w i+ i(y),h x )), where ui(w n (y),h x ) 
is defined in (25 1, ui ov (Aw,wi,W2, h) is defined in 
(26 1 and 7?.(y) is the slice of the map support TZ at 



the ordinate y. Similarly, the contribution from the i- 
th parallelogram Q ov .i to A x (h x ,lZ(y)) is given by 



F. MONTORSI et al: ON THE PERFORMANCE LIMITS OF MAP- AWARE LOCALIZATION 15 



{1,3,5} 

{2,4} 



Pi 




Figure 9. Representation of the set Q and of its subsets {Qi} igJ ^ and 
{Sou,i}jgj\/i! in the (t,u) plane (N r = 3 is assumed). Note that the 
contribution from the "discarded" parallelogram is neglected in the evaluation 
of both {24}-(23} and gl). 



u) ov (Awi(y),Wi-i(y),w i+ i(y),h x ). Therefore, the overall 
contribution to \ x (h x ,lZ(y)) can be expressed as 

} J uj ov (Aw n (y),w n -i(y),w n+ i(y),h x ) (56) 

whereas that to "f x (h x ,7Z.(y)) as 

2 ]T ^{\^n{y),h x \ (57) 



Finally, integrating (56 1 and pfy produces (23i and (24i, 
respectively. 
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